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Abstract. Seismology of stars that exhibit solar-like oscillations develops a growing interest with the wealth 
of observational results obtained with the CoRoT and Kepler space-borne missions. In this framework, 
relations between asteroseismic quantities and stellar parameters provide a unique opportunity to derive 
model-independent determinations of stellar parameters (e.g., masses and radii) for a large sample of stars. I 
review those scaling relations with particular emphasis on the underlying physical processes governing those 
relations, as well as their uncertainties. 
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1 Introduction 



The determination of accurate stellar parameters (mass, radius, effective temperature, age, and chemical com- 
position) is a fundamental and longstanding problem in astrophysics (e.g. Soderblom|2010 ). When compared to 
many other methods, seismology is recognized to provide the most precise determination of those fundamental 



parameters (e.g. Lebreton 2011). Nevertheless, such a determination is only possible by means of the use of 
stellar models and therefore suffers from our deficient knowledge of the physical processes taking place in stars 
(e.g. |Goupil et aL|2011b|a| ) 



The situation recently improved with the advent of space-borne asteroseismology, more precisely with the 
launch of CoRoT flBaglin et aL]|2006a|b| |Michel et aT]|2008| ) and Kepler ( |Borucki et~aT||2010| >. Those two 
spacecrafts are providing us with high-quality photometric data. Up to now, several hundreds of main-sequence 
stars with solar-like oscillations have been detected and several thousands oscillating red-giant stars, allowing 
for statistical analysis. With such large number of stars it is not possible to performe classical seismology, i.e. 
by individual mode fitting of the power spectrum. This is very time (and man-power) consuming so that a 
new method emerged through the use of seismic global parameters. The latter are typical global characteristics 
of the oscillation spectra such that the regularities in frequency (or period) , or the frequency of the maximum 
amplitude. 

This approach gave birth to the ensemble asteroseismology, which cornerstones are the relations between 
global seismic quantities and stellar parameters. It allows ones to infer model-independent stellar parameters 
as well as information on stellar structure and evolution. Scaling relations between asteroseismic quantities 
and stellar parameters such as stellar mass, radius, effective temperature, and luminosity have initially been 
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Brown et al. 


1991 



using ground-based data. CoRoT and Kepler confirmed these results by providing accurate and homogeneous 
measurements for a large sample of stars from main-sequence to red-giant stars (e.g., Mosser et al.|2010 ; Baudin 
et al.|2011a| |Mosser et aT|2011a|b[ |2012b|a| |Samadi et aL|2012[ ) 



Among them, the relations between the large separation (Aj/) and the mean density (p) as well as between 
the frequency of the maximum height in the power spectrum (i/ max ) and the photospheric cut-off frequency (v c ) 
have been known for a long time in the context of solar-like pulsators (e.g., Ulrich 1986 Brown et al. |1991| 



Kjeldsen & Bedding 1995b Belkacem et al.|[2011 ). In fact, as shown below, the physical grounds of these two 



relations were known for an even longer time in the context of classical pulsators (stars exhibiting opacity-driven 
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modes). This will be the subject of the first part of this review, with particular emphasize on the theoretical 
side. We will discuss the fundamental physical concept underlying these scaling relations and show how one 
can use them to infer stellar masses and radii. Section 3 will be dedicated to scaling relations that exhibit an 
important potential in providing information on the core properties of stars, or effective temperature. Last but 
not least, we will discuss the uncertainties on the seismic global parameters by showing that many biases can 
exist and are still to be addressed. 

Finally, we note that this paper is by no means an exhaustive review. Rather it puts forth the necessity of 
understanding the theoretical ground of the scaling relations, unfortunately too often bypassed in the recent 
literature. 



2 The canonical scaling relations 

In this section, we will first address the problem of the relations between the large separation (Av) and the 
mean density (p), as well as the frequency of maximum height in the power spectrum (^ max ) and the cut-off 
frequency (v c ). Those two scaling relations provide an estimate of the mass and the radius. They are now 
widely used and we will see that the underlying physics is known for a long time. Consequently, in this review 
we will denote them as canonical scaling relations. 



2.1 Relation between the large separation (Av) and the mean density (p) 



Let us first define the large separation (Av). To this end, it is necessary to introduce the first-order asymptotic 
relation that permits us to express modal frequencies (v n ,e) as a function of the structure of the star, i.e. 



R dr 



(2.1) 



where n is the radial order, I the angular degree, a a term that accounts for the near-surface effect (e.g. 
Christensen-Dalsgaard &: Thompson||1997 Rosenthal et al. 1999), and c s the sound speed. Note that such an 
asymptotic analysis assumes that we consider high radial-order modes (see Tassoul|[l980 ). 

Hence, from Eq. (2.1), the large separation which is defined as the frequency separation between two con- 



secutive radial orders (for a given €) is given by 



Av = 



Vn+\,t. — v nJ 



dr 



(2.2) 



and therefore represents the inverse of twice the time for a perturbation of pressure to cross the entire star (in 
other words a back and forth of a long wavelength pressure wave) . 

Now, we aim at deriving the scaling relation between Av and the mean density. To this end, one has to 
adopt the homology relations. So let us consider two stars such that for two shells verifying r/R = r' /R', the 
corresponding mass shells equal (m/M = m'/M'), where M,M' are the total masses of two stars belonging of 
a homologous series, and R, R' their total radii. This type of approximated stellar models has been extensively 
discussed in the literature (e.g. Cox &; Giuli|1968| |Kippenhahn &; Weigert|1990 ), and it is possible to show that 
pressure (p) and density (p) of both models are related by (see Sect. 20.1 of Kippenhahn fc Weigert|[l990 ) 
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Then, from Eqs. (2.3) and (2.4), the relation of the sound speed immediately follows 

-1/2 
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(2.3) 
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(2.5) 



Now, to demonstrate the relation between the large separation and the mean density, let us define the ratio 

Av 



K 



Av 1 



(2.6) 
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Using Eq. (2.5) together with the relation r/R= r'/R', it is straightforward to demonstrate the desired scaling 
relation, i.e. 

-i -l 
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R dr 
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(2.7) 



In other words, the large separation of a given star (Ai/) can be related to its mean density such as 

1/2 



Av 



^ A, , 
P&J 



(2. 



where, as classically found in the literature, the Sun has been used as the reference. Equation (2.8) demonstrates 
the scaling relation between Av and p, and shows that the only underlying hypothesis is the homology, which 
will be discussed in Sect. 14.2.11 

At this stage, it is worthwhile to recall that Eq. (2.8) and its derivation is not a novelty and was known, in 
a hardly different framework, for many decades. Indeed, for classical pulsators, the mode frequencies are often 
near the fundamental mode frequency (e.g., Cepheid, (5-Scuti, /3-Cephei). For fundamental radial mode (vq), its 
period (IIq) is therefore proportional to the time for a pressure perturbation to cross the entire star, i.e. 



n oc 



dr 



n oc p 



1/2 



(2.9) 



This relation (Eq. 2.9 1 is quite famous in the context of classical pulsators and one of the first authors to mention 



it was, to our knowledge, |Eddington| ( |1917[ ) for explaining the periodic motion of Cepheids. Subsequently, the 
relation between 
stars as pulsations 



relation between vq and the mean density had been often used as an argument to identify periodic motions of 
Moreover, it is the basis of the famous period-luminosity relation of Cepheids (e.g 



et al. 


1972 


Cox 


1980) 



and was proposed for a long time by several authors such as Ledoux & Walraven ( 1958 ) (see Cox 1980 for a 
comprehensive discussion) . 



As a conclusion, the derivation of Eq. (|2.8|) is nothing but the same as for Eq. (2.9) 



Cox 



2.2 Relation between the frequency of the maximum height in the power spectrum (V max J and the cut-off frequency 

fa) 

The derivation of the relation between v max and v c is more difficult. Because sometimes not fully understood 

), we explain here how the basic physical picture is grasped and that 
departure from the observed relation arises from the complexity of non-adiabatic processes involving time- 
dependent treatment of convection and not from the failure of the physical picture. In addition, a discussion on 
the physics of opacity-driven pulsations in stars will bring us with the conclusion that the main physical reason 
for the existence of the relation between ^ max and v c is a common feature of pulsating stars. 



(e.g., Christensen-Dalsgaard 2011a|b 



2.2.1 Derivation of the v„ 



relation 



Let us first begin by recalling that solar-like oscillations are the result from a balance between mode driving 
and damping. Therefore, as a first approximation, each mode can be considered as a driven and damped 



oscillator. The driving is related to turbulent Reynolds stresses (see Samadi 2011 for a comprehensive review 
on stochastically excited modes) while damping is caused by a combination of physical processes not discussed 
here (but see |Houdek| < p008| ) for a discussion on mode damping). 

Therefore, the frequency ^ ma x is in principle determined by both physical mechanisms. To have a more 
precise view of what governs i/ max , one has first to determine which of the aforementioned process is responsible 
for the maximum height in the power spectrum, thus the physical mechanism controlling f max - We thus consider 
the height H of a given mode in the power spectrum, which is a natural observable. For stochastically excited 



modes, the height of the mode profile in the power spectrum is (e.g. Baudin et al. 2005 Chaplin et al. 2005 
Belkacem et al.||2006| ) 

H = Wm > (2 ' 10) 
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where P are the excitation rates, r\ the damping rates, and Ai the mode masse^] It turns out that the depression 
(plateau) of the damping rates 77 is responsible for the presence of a maximum in the power spectrum. This 
is in agreement with theoretical computations (Houdek et al. 1999 Chaplin et al.|[2008 Bcl kacem et al.pOll 



2012| ) and recent observations of the solar-like stars by Kepler ( Appourchaux et al.||2012 ) 



The subsequent issue then relies on the physical origin of the depression of the solar damping rates. This 
plateau originates from a destabilizing effect in the super-adiabatic layers and occurs when the modal period 
nearly equals the thermal time-scale (or thermal adjustment time-scale) in the superadiabatic layers. This was 



first mentioned by Balmforth (1992) (see his Sect. 7.2 and 7.3) and confirmed by Belkacem et al. (2011). The 



authors used two different non-adiabatic pulsation codes, making this conclusion quite secure. This can be 
expressed by the following condition 



Q = 2wv n 

with the inverse of the thermal time-scale defined as 



1 



L 

Airr 2 pc v TH p 



=T~ 1 +T- 1 
conv 1 rad > 



(2.11) 



(2.12) 



where L is the luminosity, r the radius, p the local density, c v — (dU/dT) p with U the specific internal energy, 
Hp the pressure scale height, r ra( j and T con v the radiative and convective thermal time-scales, respectively. Note 



that Eq. (2.12) is a local formulation of the thermal time-scale; a non-local one can also be defined (e.g. Pesnell 



1983). It is important to note that, contrary to the situation for classical pulsators, for which the envelope 



is dominated by the radiative transport of energy, Eq. (2.12) exhibits contributions of both the radiative and 
convective fluxes. 

The last step is to establish the relation between the thermal time-scale (r) and the cut-off frequency. In 
the mixing- length theory framework the thermal time-scale can be recast such as ( Belkacem et al.||2011 ) 



conv 

t pc v TH p 



ad 



Mi 

a 



2Hr, 



mnv 



(2.13) 



where M. a = v conv /c s is the Mach number, and a the mixing-length, Fconv>-frad the convective and radiative 



fluxes respectively, c s /(2H p ) the cut-off frequency. Note that the cut-off frequency appears in Eq. (2.13) as the 
ratio c s /(2H p ). 

This relation is verified, in Fig. [T] by using a grid of stellar models. The relations between the thermal 
time-scale and the cut-off frequency is very tight for red giants while a dispersion is observed for main-sequence 
stars. Such a dispersion is explained by the dependence to the Mach number to the third power. Indeed, this can 
be easily understood since the Mach number predominantly depends on the effective temperature that varies 



much more for main-sequence and sub-giants than for red giants. Equation (2.13) is useful to explicitly show 



the relation between r and v c , but an investigation of this relation using a set of 3D hydrodynamic numerical 
simulation would be desirable in the future to get more quantitative estimates. 

It is nevertheless worth to emphasize that the resonance betw een th e local thermal time-scale and the 
pulsation period (which is nothing but the transition region, see Sect. 2.2.2) explains the fc^/i^ observed relation 
between ^ max and v c . More interestingly, as observed in Fig. [2j the maximum dispersion predicted for main- 
sequence stars is in agreement with recent observations of Bedding (2011). Indeed, as depicted by Fig. [2j 
a departure from the scaling relation seems to occur for main-sequence stars (i.e., large i/ max ) in qualitative 
agreement with Fig. [T] 



2.2.2 The transition region: a common feature of both opacity-driven and solar- like pulsators 

As explained in the preceding section, the cornerstone of the physical picture underlying the v ma ^-v c relation 
is the nearly equality between the thermal time-scale and the modal period in the super-adiabatic layers. 



Let us first briefly introduce the notion of transition region (see Cox 1974 1980, for a more complete 
discussion). Such a region is defined by the layer for which the ratio of the thermal time-scale to the modal 
period equals unity. Indeed, it separates the quasi-adiabatic from the non-adiabatic layers, more precisely 



*The mode mass corresponds to the total amount of mass effectively moved by a given mode. Its definition is M = 
f™\Z\ 2 drn/\t(r = R)\ 2 . 

T Note that the notion of tightness sometimes found in the literature is very subjective and would need a more rigorous quantitative 
estimate. 
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Fig. 1. Thermal frequency (1/r) computed from Eq. ( 2.13 1 versus the cut-off frequency (computed as the ratio c s /(2Hp)), 
normalized to the solar values, for models with masses ranging from M = 1.0 Mo, M = 1.2 Mq and M — 1.4 Af@ from 
the ZAMS to the ascending vertical branch. The inputs physics of the models can be found in Belkacem (2011). 
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Fig. 2. Observed ^ max versus predicted i/ max , computed using Eq. (2.171, for stars observed from the ground (see Bedding 
[2011] for details). Figure from |Bedding| poTT| . 



• In quasi-adiabatic regions, there is a negligibl^] exchange of energy between the oscillations and the 

tBy negligible we mean on the mode stability. 
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background. In those regions, Q^> 1 (see Eq. 2.11 1 then r » II osc , so that the background does not have 
time enough to react to the perturbations introduced by one cycle of the oscillation. Consequently, over a 
time-scale r perturbations related to the oscillation compensate and the background experiences a nearly 
vanishing average over many oscillation cycles. 

• In non-adiabatic regions, the energy exchange between the oscillation and the background dominates. 
We are in the opposite situation where Q <C 1 so r <C n osc . Indeed, in this case the thermal structure 
"instantaneously" adjusts to the perturbations generated by the oscillation so that, say during compression, 
any accumulated energy is immediately lost and the perturbation of luminosity is nearly constant. 

This transition region is an essential ingredient of the K-mechanism in opacity-driven pulsators. In those stars, 
pulsations are destabilized by the perturbation of the opacity. But to be efficient, this destabilization must 



fulfill several conditions (e.g., Cox 1980 Cox & Giuli 1968 Pamyatnykh 1999) among which the transition 



region must lie in ionization region. Indeed, the destabilization occurs in the ionization region and if r <C n osc 
the thermal structure adapts so quickly that the flux is frozen. In the opposite situation if t 3> n osc we are in 
the quasi-adiabatic situation as described above. Finally, one must have r ss II osc in the ionization region for 
the destabilization to be efficient and dominant over damping terms. 

The situation is similar in solar-like pulsators, except that the destabilization by the perturbation of the 
opacity never dominates over damping terms and the situation is complicated by the presence of convection 



which modifies the thermal time-scale (see Eq. 2.11|. As shown by Balmforth (1992) and Belkacem et al. (2011) 



the depression of the damping rates occurs when the destabilization by the perturbation of the opacity becomes 
efficient (and partially compensates the other damping mechanisms), i.e. when the transition region nearly 
coincides with the ionization region^] Indeed, for low-frequency modes, Q > 1 so that the destabilization is 
inefficient and, for high-frequency modes, Q> 1 the important destabilization in the super-adiabatic layers are 
compensated by the damping in the atmospheric layers. Therefore destabilization have its maximum impact 
on the total damping rate for Q ss 1 in the super-adiabatic layers (near the ionization region of hydrogen), in 
other words when the transition region occurs in the super-adiabatic layers. 

We then conclude that the observed i/ max - v c relation simply results from a common feature of pulsating 
stars, i.e. the occurrence of the non-adiabatic effects in the transition region. 



2.3 Inferring seismic masses and radii from global seismic quantities 



Given the scaling relation described in Sects. |2.1| and |2.2| it is now possible to link seismic global quantities to 
stellar parameters. 

Let us start with the scaling relation relating the large separation to the mean density of the star (Eq. [2? 

i.e. 



Ap oc p 1 / 2 oc 



M 



1/2 



(2.14) 



where Av is the large separation, p the mean density, M the total mass of the star, and R its total radius. 
The second one relates the frequency of the maximum height in the power spectrum to the cut-off frequency, 

i.e. 

(2.15) 



The cut-off frequency can be approximatively defined for an isothermal atmosphere (see Stello et al. 2009 for 
a discussion) such as 

v c = — , (2.16) 

where c s is the sound speed. Through the hydrostatic equilibrium, the pressure scale height is related to the 
gas pressure by P = pgH p and using the ideal gaz equation of state P oc pT, Eq. (2.16) becomes 



v c (X 



9 M 

Vt r 2 Vt 



(2.17) 



§Note that in solar- like pulsators, the hydrogen ionization region is located in the uppermost atmosphere and nearly coincides 
with the super-adiabatic region 
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From Sect. |2.2| it becomes clear that the introduction of the cut-off frequency is more historical than physically 
justified. Indeed, Brown et al. ( 1991 ) conjectured a relation between f max and v c but from a physical point of 



view it would be more rigorous to mention the relation between v max and the thermal time-scale (Eq. 2.11), 



even if the cut-off frequency can be artificially introduced as in Eq. (2.13). 

Now, using Eqs. (2.14) and (2.17), one obtains the desired relations (normalized to the solar values) 



M 



oc 



max 



R 



max 



Av 
A^ 



T e ff,Q 
eff,© 



T. 



3/2 



1/2 



(2.18) 
(2.19) 



where T c g is the effective temperature. Note, however, that from Eq. (2.17) to Eqs. (2.18) and (2.19) one 



implicitly assumes that T = T e g. This is not obvious and would deserves further investigations. 



Equations (2.18) and (2.19) constitutes the main frame of what is now commonly called the ensemble 



asteroseismology and recently gives rise to an important work that provide us the tools for a new grip on stellar 
physics. An all-comprehensive review of the work that make use of the above-mentioned scaling relations would 
be tedious but it is worth to emphasize the diversity of the applications, namely 

• Model- independent determination of stellar parameters. Such a type of application is now commonly 
used to infer masses and radii, so as to identify the observed stars (e.g., |Kallinger et al. 2010 Mosser 



et al. 20101, but also for other applications such as the characterization of planets (e.g. 
2012| Jackiewicz et al.||2012[ ). 



Borucki et al. 



• Constraint on stellar evolution and population. From the determination of stellar parameters for a large 
set of stars, ranging from the main-sequence to the red-giant branch, it is obvious that one of the first 
applications is to give constraint on stellar evolution and populations. Such work has been recently 



performed by many authors among which Mosser et al. (2012a I for constraining mass loss at the tip of the 



(2012) for constraining populations in the milky 



red giant branch, or |Miglio et aT] ( |2009| |2012a[ ); |Miglio 
way. 

A distance indicator. From the knowledge of stellar radius and with the effective temperature as an input, 
the distance is derived from a comparison with apparent magnitudes. It is claimed that such a method 
can provides accurate results, of the order of 10% ( Miglio et al.||2012b Silva Aguirre et~aI1[20T2| . 



• Improved determination of log g and T e g. A striking example of the use of seismic scaling relations is the 
determination of stellar surface gravities. Classically, surface gravities are obtained by using spectroscopy 
and isochrone fitting. However, a look at Eq. (2.17) immediately shows that the relation between ^ max and 
v c gives access to the surface gravity. It has been shown bylMorel & Miglio (2012 1 that the seismic gravity 



is compatible with classical methods and that one can expect a much better accuracy from seismology than 
from the other methods. Last but not least, the use of this seismic gravity as an input in spectroscopic 



analysis provides better determination of effective temperature (e.g. Batalha et al. 2011 Creevey et al. 



2012) 



Those few examples offer an overview of the current multiple use of the scaling relations, but ensemble 
asteroseismology is not limited to these few examples and other scaling relations are promising. 



3 Some promising scaling relations 

In this Section, we will address three relations that from our point of view are very promising in providing us 
highly valuable information for stellar physics, namely the age, luminosity, and effective temperature of stars. 
Those relations are just beginning to be exploited, for different reasons we will discuss below. 



"It would be more accurate to replace the term model-independent by stellar-model-independent since the derivation of the 
scaling relations rely on physical assumptions, therefore on a modeling. 
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Fig. 3. Gravity- mode period spacing as a function of the pressure-mode large frequency spacing. RGB stars are indicated 
by triangles; clump stars by diamonds; secondary clump stars by squares. Small gray crosses indicate the bumped periods 
measured by Mosser et al. (2011a). The solid colored lines correspond to a grid of stellar models (see Mosser et al.|2012b 
for details). Figure from Mosser et al. ( 2012b I. 



3.1 The All - evolutionary status relation 



Up to recently, mainly pressure modes have been detected. In contrast, gravity modes (whose restoring force 
is dominated by buoyancy) were very difficult to detect, especially in the solar case since they have very low 
amplitudes (e.g., |Belkacem et al.||2009{ Appourchaux et aL||2010 ). Nevertheless, the situation changed with the 
advent of space-borne missions. Indeed, since the first unambiguous detections and identification of non-radial 
modes in red giants by CoRoT ( |De Ridder et al.||2009 ), a great leap forward is being experienced by the stellar 
physics community. This has been possible because of the large amplitudes that oscillations develops in red 



giants compared with the main-sequence stars (Baudin et al. 2011a Samadi et al. 2012). It then makes the 



detection easier, allowing for the identification and characterization of a very large number of red giants (e.g., 
Mosser et aT|2010 |. 

In main-sequence stars, p-modes have high frequencies while <?-modes have low frequencies. In contrast, in 
red-giant stars both p and g modes lie in the same frequency range. This is due to the structure of red giants, 
since the radiative core of red giants contracts from the end of the hydrogen burning phase while the envelope 
expands. Therefore, the total radius increases, the mean density decreases, then p-mode frequencies decrease 
during the evolution (see Eq. (2.9)). For g modes, we are in the opposite situation since their frequencies 
increase as the result of the contracting core and more precisely due to the increase of the buoyancy frequency 
(see Eq. (3.1 )). Then, on the red giant branch, the p and g modes are in the same frequency range, then modes 



propagate in both the outer and inner cavities. Moreover, these cavities are coupled by an intermediate zone 



in which modes are evanescent. Consequently, red giants exhibit what we call mixed modes (e.g., Dziembowski 
et al. 2001 Dupret et al. 2009), i.e. with a g nature in the core and a p nature in the envelope. As a result, 



while probing the core they have enough amplitude at the surface of the stars to be detected. Those peculiar 
physical properties of the oscillations of red giants lead to the detection and characterization of mixed modes 
in a large sample of stars (e.g., Bedding et al.|20lT Mosser et al.||2011a 2012b), giving us the opportunity of a 
new grip on stellar physics of advanced evolutionary stages. 

To go further, it is first essential to introduce the period spacing. It is the counterpart of the large separation 
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for gravity modes. Indeed, in the asymptotic regime, gravity mode periods follow the relation ( Tassoul]|1980 ) 

n„, f = An ( (n + £ 9 ), (3.1) 



where e g is a phase shift, and AIT the period spacing given by 



An, = 



2tt 2 



ri 



r 



(3.2) 



where r^ and r 2 are the radius of the inner and outer turning points of the g-mode cavity. It is worth to note 
that the integral in Eq. (3.2) is related to the evolutionary state of the star. Indeed, as the star quits the 



main-sequence, its core contracts and the buoyancy frequency N increases, leading to a decrease of the period 
spacing. Therefore, the detection of mixed modes and the period spacing from the observation permitted to 



assess the evolutionary stage of red-giant stars as shown by Bedding et al. (2011) and Mosser et al. (2011a 



More precisely, low-mass stars on the ascending red-giant branch (RGB) and after helium ignition (clump stars) 
can be found at the same location in the HR diagram. Since such stars have similar envelopes it was impossible 
to differentiate them from the use of p modes only. In contrast, the detection of mixed modes and the period 
spacing makes it possible to distinguish between those two evolutionary stages as shown by Fig. [3j 

Nevertheless, there is a subtlety that cannot be avoided. In fact, the period spacing between two observed 
mixed modes does not follow exactly the asymptotic relation given by Eq. (3.2). This is related to the mixed 



nature of the mode, since the observed modes are not pure <?-modes they are affected by their acoustic nature 



and a departure from Eq. (3.2) follows. To cope with this issue, Goupil (20121; Mosser et al. (2012b), proposed 
an asymptotic relation adapted for mixed modes (based on the formalism developed by Shibahashi (1979])), i.e. 



Au 

-i H arctan 

7T 



gtan7r 



/ 1 



\An, 



(3.3) 



where v m is the mixed mode frequency, f np ,£=i the frequency of pure p modes, and q a coupling factor. Therefore, 
using Eq. ( |3.3[ ) with the measured v m yields a determination of the period spacing (AIT). This is a crucial step 
for comparing the observations and the modeling as well as for determining the evolutionary stage, as illustrated 
by Fig. 

Finally, one can conclude that the scaling relation between the period spacing and the evolutionary state of 
stars is only in its infancy but still very promising since it provides unprecedented information of the innermost 
layers of stars. 



3.2 Scaling relations related to mode amplitudes and linewidths 

Asteroseismology is not limited to the analysis of mode frequencies but also to mode amplitudes and linewidths. 
The latter being related to the exchange of energy between the oscillation and the background, and not only 
to the star structure, they are subject to the uncertainties related to the coupling between convection and 
pulsation. Nevertheless, several successful attempts have been proposed to scale those seismic parameters to 
the stellar fundamental parameters. 



3.2.1 Mode amplitude vs L/M 

We first consider the relation between mode amplitudes and stellar parameters, both in terms of velocity and 



intensity fluctuations. On the basis of the theoretical calculations of Christensen-Dalsgaard & Frandsen (1983), 
Kjeldsen & Bedding ( 1995a| ) derived the first example of a scaling relation given in terms of the maximum of 
the mode surface velocity (hereafter Vmax). This scaling predicts that 



(3.4) 



where L is the luminosity. Christensen-Dalsgaard & Frandsen ( 1983 1 assumed that there is an equipartition 
between the energy carried by the most energetic eddies and the modes. As mentioned by |Belkacem et al.] 
( |2009[ ) and |Samadi| ( |2011[ ), a necessary (but not sufficient) condition for having such an equipartition is that 
turbulent viscosity is the dominant source of damping. However, there is currently no consensus as to what is 
the dominant physical processes contributing to the damping of p-modes. 
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Prior to the CoRoT mission, observations of mode velocity in solar-like stars were sparse (e.g. Houdek 
& Gough||2002|) but motivated several theoretical work on the physical mechanisms underlying mode driving 



( |Houdek et al.|1999| |Houdek fc Gough|2002| |Samadi et al.|2007|[Samadi|2011| ). With the launch of space-borne 
mission CoRoT, and its unprecedented high-quality data, such work on scaling relation of mode amplitudes 
became achievable (Michel et al. 2008; Bau din et al.|2011a |. This has been confirmed with the Kepler mission, 
which motivated an important work on those relations (e.g. Kjeldsen & Bedding 2011 Huber et al.|2011 Mosser 
et al.||2012a[ ). 

The large amount of high-quality data from CoRoT and Kepler led to a variety of results that paradoxically 
complicated the picture since no clear scaling relation emerged. Samadi et al. (2012) addressed this issue for 



red giants, based on theoretical developments, a set of 3D numerical simulations, and the observations. They 
found that the maximum amplitude (in term of velocity) follows the scaling relation 



Vr 



max '/max 



1.3 



M 



0.525 



(3.5) 



where rj max is the linewidth at v — v max . In terms of intensity fluctuations, Samadi et al. (20121 showed 
the necessity to go beyond the adiabatic relation between velocity and intensity, especially for red giants, and 
proposed 

Y \ 1.55 / i\ r \ 0.5 

' ' ' 1 (3-6) 
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These scaling relations show a systematic discrepancy for red-giants, which is attributed to non-adiabatic effects. 
Consequently, we can state that the scaling relation of mode amplitudes is now well understood for main-sequence 
stars but still need to be investigated for red giants. 



3.2.2 Mode linewidths vs T cS 

For mode linewidths (or equivalently mode damping rates), scaling relations have been investigated only very 
recently. This is the result of the need for long-time and almost-uninterrupted monitoring to resolve individual 
modes and to enable their precise measurements. 

Houdek et al. (1999), and later Chaplin et al. (2009), have investigated the dependence of mode-damping 



rates on global stellar parameters. From ground-based measurements, Chaplin et al. (2009) found that observed 
mode linewidths follow a power-law of the form rj cx T^ s (where T e s is the effective temperature) and no clear 
tendency emerged when rj is scaled with the ratio L/M. Nevertheless, these measurements were based on 
short-term observations and derived from an inhomogeneous set of analysis and instruments, resulting in a large 
dispersion. This was settled by Baudin et al. ( 2011a|b ) (Fig. [3]) using a homogeneous sample of CoRoT data. 
They found that a unique power-law hardly describes the entire range of effective temperature covered by main- 
sequence and red-giant stars and proposed that mode linewidths of main-sequence stars follow a power-law of 
rjg , while red-giant stars only slightly depend on effective temperature (T G ^°' 3±0 ' 9 ). The latter result was later 
confirmed and extended by Kepler observations (Fig. [4]) to main-sequence and sub-giant stars ( Appourchaux 
et al.|2012 ). We also note that Corsaro et al. (2012) proposed that mode linewidths follow an exponential power 
law. In absence of a strong theoretical argument to adopt either a power law or an exponential, the statistical 



significance must dictate our choice and this is still to be performed. From a theoretical point of view, Chaplin 



et al. (20091, based on the formalism developed by Balmforth (19921; Houdek et al. (1999) and Chaplin et al. 



(2005), were ab 



(2005 ), predicted a power-law of r\ tx T e 4 ff which disagrees with CoRoT and Kepler observations ( Houdek|2012 ). 



e to reproduce 



In contrast, Belkacem et al. (2012), based on the formalism of Grigahcene et al. 
both CoRoT and Kepler observations. 

To get more insight into the relation between r\ and T e g, let us first write down the integral expression of 
the damping rates (e.g. Grigahcene et al.||2005 ) 
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(3.7) 



where uj is the mode frequency, 5p the Lagrangian perturbation of density, SP the perturbation of the total 
pressure (including the turbulent pressure), po the mean density, and where the star denotes the complex 
conjugate. The mode inertia is given by 

(3.8) 
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Fig. 4. Mode linewidths (normalised by the solar value, r sun = 0.95 /iHz) versus effective temperature. The squared 
symbols represent theoretical calculations computed as explained in Belkacem et al.| ( |2012[ ) . The triangles correspond to 
the observations of main-sequence stars derived by Appourchaux et al. ( 2012[ ) from the Kepler data (with their 3-<r error- 
bars). The dots correspond to the observations of red giants (with T e ff < 5200 K) and main-sequence (with T e a > 5200 
K, with their 3-<r error-bars) stars as derived by Baudin et al. ( 2011a|b I from the CoRoT data. 



where £ is the eigendisplacement vector. Therefore, a look at Eqs. (3.7 1 and (3.8) makes it clear the need for 
disentangling the effects of mode inertia and the work integral (i.e., the integral appearing in Eq. |3.7[ ). The latter 
is related to non-adiabatic processes corresponding to a transfer of energy between pulsation and convection. 
Hence, it can be assumed at first glance that the work integral scales dimensionally with the ratio L/A7|^[ As 
verified by Belkacem et al. (2012), it follows that the relation 



7] I oc 



2.7 



(3.9) 



holds. In contrast, the mode inertia (J) does not depend on mode energy leakage but on the star's static 
structurj**] and more precisely on the properties of its uppermost layers . Henc e, one can expect mode inertia 



to scale with the surface gravitjEM It has been shown in Belkacem et al. (2012) that 



I oc g 



-2.4 



(3.10) 



Using Eq. (3.9) and Eq. (3.101, it turns out that 



This simple analysis allows us to explain qualitatively the strong dependence of mode damping rates on effective 



temperature. Finally, it is important to stress that Eq. (3.11) is compatible with current K epler observations 
(T. Appourchaux, private communication). 



4 Do seismic scaling relations give us accurate stellar global parameters? 

This is a difficult and still open question issue. We explain here what are the main sources of uncertainties. 
The latter includes biases, all unknown and missing processes, as well as the precision of the measurements. 



II This comes from the perturbed energy equation that shows that entropy perturbation dimensionally scales as the ratio L/M. 
"Except through non-adiabatic effects on mode eigendisplacement. 

Note that mode inertia also scales with the dynamical timescale ^J~{GM/B?) with almost the same dispersion as for the surface 
gravity. 
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The typical precision of the stellar seismic indices is very good (better than 1% for Av and 5% for v max in 
most cases). Under the assumption that the effective temperature is determined with a precision of, say, 100 
K, it translates into a precision of 20% for the mass, 8% for the radius, and 0.04 dex for logg. It is worth 
to mention that the precision on v max is the limiting factor. These numbers are very encouraging and often 
presented as the uncertainty, but they do not include biases. Indeed, efforts are currently undertaken to assess 
the reliability of those numbers and more precisely to identify the biases from both the observational and the 
theoretical side. 



4.1 Observational uncertainties 

To illustrate the importance of considering the possible sources of biases with great care let us consider two 
examples, namely Av and z/ max . 



4.1.1 Observational determination of Av 



As demonstrated by Verner et al. (2011 ) on the basis of hare-and- hound exercises and fitting of Kepler data by 



six teams, the way the large separation is determined plays an important role. Indeed, the authors conclude 
that the expected relative precision on Av is about 2%. 

We also emphasize that the large separation is a quantity that is relevant in an asymptotic regime only. Thus, 
it can be extracted if the assumption underlying the asymptotic analysis is respected, here the requirement is 
that one must consider high radial-orders. However, in practice, solar-like oscillations are observed for low to 
moderate values of the radial orders, especially for red giants (n £ [6; 15]), which make the basic assumption 
violated. This is a non-negligible source of biases as demonstrated by Mosser et al. (2012c I 
systematics as high as of 8% for the mass and 4% for the radius determination. 



and can generate 



4.1.2 Observational determination of v„ 



The second example concerns the determination of v max , for which several issues arise. Since we are dealing 
with a stochastic process, the maximum of height in the power spectrum can vary depending on the considered 
time-duration of the observations. This makes the determination of y max quite unsecured and requires very long 
observations to settle the problem (unfortunately this is not always possible) and to approach the stationary 
state. Therefore, except if we are able to demonstrate that stochastic effects are negligible, it is safe to consider 
that a possible bias of about half a large separation. The related uncertainties can be roughly estimated by 
considering the ratio Av/(2v mayi ). It is about 2.5% for main-sequence stars and 5% for red-giants. Note that 
such an effect has been considered for red giants by Mosser et al. (2011b). In addition, as for Av, several 
methods exists for deriving ^ max from the light-curve and Verner et al. (2011) found an average precision of 
about 4%. 

Moreover, another issue naturally arises; can v max be observationally determined by both considering the 



maximum of amplitude or height in the power spectrum? Following the work of Belkacem et al. (20111 the 



suitable choice (in the sense it is physically grounded) is to use the maximum height. Indeed, the relation 
between v max and v c is due to the occurrence of the transition region in the superadiabtic layers that translates, 
from an observational point of view, by the depression of the damping rates. Hence, the quantity one must 
consider must be dominated by the damping processes. It is the case for mode height, but not for mode 
amplitude. Indeed, mode amplitude derives from a more subtle balance between mode driving, damping, and 
mode inertia. Therefore, the maximum of height and amplitude is expected to be different. 



4.2 Theoretical uncertainties and additional dependencies 

As illustrated in Sect. |4.1[ biases can arise from the observational determination of seismic indices. However, 
as shown in Sects. |2.1| and |2.2[ scaling relations are based on modeling and therefore cannot be considered as 
perfect, so does the derived stellar parameters. 

In the following, we discuss briefly the main physical reason underlying the dispersion of those relations. 



4.2.1 The limits of the homology assumption and the Av - p relation 

As depicted in Sect. |2.1[ the physical assumption that permits us to relate the large separation to the mean 
density of a star is homology. While often considered as a crude approximation to derive the internal structure of 
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Fig. 5. Ratio between the large separation and the squared mean density of stars, normalized to the solar values, from 
the ZAMS to the ascending red-giant branch (see White et al.|2011 for details). Figure from White et al. (2011 1. 



a star, this approximation nevertheless gives ones quite a good result concerning this scaling. Figure [5] illustrates 
that; for typical solar-like oscillating stars, this relations holds within 3 to 4%. 

To understand such a departure, it is necessary to recall the main physical hypothesis that make the homology 
strictly valid. The first main requirement is hydrostatic equilibrium, i.e. that the acceleration term in the 
momentum equation must vanish. The second one is thermal equilibrium, i.e. that the energy generated must 
be strictly compensated by the energy emitted (in other words, dS/dt = 0, where S is the entropy). Finally, as 



demonstrated by Cox & Giuli (1968), the homology requires that the constitutive equations (such as opacity, 
production rate of energy, etc..) must be power-laws of their arguments. 

It is clear that during its evolution, a star breaks all the requirements implying a departure from the 
homology. This scaling relation would then need a deep investigation of the physical reason explaining precisely 
the origin of the departure from the homology. This requires to consider each evolutionary state (main-sequence, 
sub-giant, and red-giant phase) separately since the physics differs from one to another. Such a work is highly 
desirable in the future, to understand why this relation is so precise and to propose improvements. 



4.2.2 Influence of the Mach number to the v n 



relation 



An other example is the relation between v max and i/ c , which is not exact as shown in Sect. 2.2 More precisely, 
the relation between the thermal frequency and the cut-off frequency is not direct and many other physical 
quantities appear, among which the Mach number is the more important. As already explained in Sect. |2.2[ 
the cut-off frequency is artificially introduced for historical reasons and this subsequently leads us to introduce 
the Mach number. We recall that this number is the ratio of the turbulent (or convective) velocity to the sound 
speed and permits us to measure the degree of turbulence of a flow. 

In solar-like stars, this number is almost negligible in the interior (M. a <C 1) but increases significantly near 
the photosphere. In these layers, convection becomes inefficient and convective velocities increase rapidly over a 
relatively small radial scale to sustain the convective flux. As a result, in this region the Mach number reaches 
a maximum (which is of the order of 0.3 for the Sun). From one star to another, this number varies typically 



from 0.3 to 0.7, mainly depending on the evolutionary status and the mass of the considered stars (Houdek 
et al.||1999| ). 
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Therefore, from Eq. p7l3] >, 
is to be investigated. 



it becomes clear that such an extra-term in the relation between 



and 



Simple theory of convection such as mixing-length theories, aside from giving very 
different results for the convective velocities ( Samadi et al.|2006 l, provides us unrealistic estimates of turbulent 
velocities (e.g. |Samadi et al. 2003). Therefore, the only way to overcome this difficult problem is to use 3D 
hydrodynamical simulations to get more insight into the evolution of the Mach number across the HR diagram 
(work in progress). 



5 Concluding remarks 

In this review, we have discussed the now commonly used scaling relations (i/ max -i/ c and Ais-p) for deriving 
seismic masses and radii. We have shown that their physical justifications were already known, for a long time, 
in the context of classical pulsators. Therefore, it is striking to note that those relations derive from common 
features of pulsating stars and not solar-like pulsators only. We also discussed how the seismic masses and radii 
are derived and emphasized the numerous and important applications for improving our knowledge of stellar 
structure and evolution. We did not limited our discussion to these scaling relations but also addressed other 
scaling relations we consider as promising, namely; the relations between mode amplitudes and luminosity, 
between mode linewidths and effective temperature, and between period spacing and evolutionary status of 
stars. 

We then focus our discussion on the uncertainties related to those scaling relations. It turns out that 
uncertainties arise from both the observational and theoretical sides. Consequently, the stellar parameters as 
derived from the scaling relations also suffers from those uncertainties. Therefore, one of the crucial step to 
obtain precise and accurate seismic stellar parameters is the calibration. Preliminary work in the direction has 
been achieved (e.g. Silva Aguirre et aL]|2012 Huber et al.||2012), and give encouraging results. 



I am grateful to M.J. Goupil, B. Mosser, and R. Samadi for reading the manuscript and for many fruitful discussions. I also thank 
T. Bedding and T. White for providing Figs. [2] and [5] 
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